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ABSTRACT 

We used extensive ground-based multisite and archival spectroscopy to derive obser- 
vational constraints for a seismic modelling of the magnetic (3 Cep star V2052 Ophi- 
uchi. The line-profile variability is dominated by a radial mode (fi ~7. 14846 d _1 ) and 
by rotational modulation (P rot =3. 638833 d). Two non-radial low-amplitude modes 
(/a=7. 75603 d^ 1 and /3=6. 82308 d _1 ) are also detected. The four periodicities that 
we found are the same as t he on es discovered from a companion multisite photometric 
campaign (Handler et al.| |2012| ) and known in the literature. Using the photometric 
constraints on the degrees £ of the pulsation modes, we show that both /2 and are 
prograde modes with (£, m)=(4, 2) or (4, 3). These results allowed us to deduce ranges 
for the mass (M e [8.2,9.6] M Q ) and central hydrogen abundance (X c <G [0.25,0.32]) 
of V2052 Oph, to identify the radial orders ni=l, ri2=— 3 and ris=— 2, and to derive 
an equatorial rotation velocity v cq g [71,75] km s . The model parameters are in 
full agreement with the effective temperature and surface gravity deduced from spec- 
troscopy. Only models with no or mild core overshooting (a ov G [0, 0.15] local pressure 
scale heights) can account for the observed properties. Such a low overshooting is op- 
posite to our previous modelling results for the non-magnetic (3 Cep star 6 Oph having 
very similar parameters, except for a slower surface rotation rate. We discuss whether 
this result can be explained by the presence of a magnetic field in V2052 Oph that 
inhibits mixing in its interior. 

Key words: stars: early-type - stars: individual: V2052 Oph - stars: interiors - stars: 
magnetic fields - stars: spots 



1 INTRODUCTION 

Stellar modelling of B-type main-sequence stars is benefiting 
greatly from seismic data of /3 Cep stars. Among the vari- 
ous classes of B-type pulsators (e.g. Chapter 2 of Aerts et 
al.ll2010Tl this is the only one with members having clearly 
identified values for the wavenumbers (I, m, n) for several of 
the detected oscillation modes, which is a prerequisite for 
successful seismic modelling. 

Models based on standard input physics could satis- 
factorily expl ain se ismic data of t he P C ep stars V836 Cen 
(Aerts et al.ll2003l ; Dupret et al.| |2004|) P CMa (Mazum- 
dar et al.l l2006t). 6 Ceti (Aerts et al.l l2006f ), and 0Oph (Bri- 
quet et al.l l2007h . However, models with the same input 
physics failed to explain the much richer seismic constraints 
assembled from huge multi-technique multisite campaigns 
organised during many months for the P Cep stars v Eri 
(De Ri dder et al.| I2004J : Pa myatny kh, Handler & Dziem- 
bo wskill2004l : Ausseloos et al.l l2004T ) an d 12 L ac (Handler et 



al.ll2006l : Dziembowski & Pamvatnvkhl l2008l : Desmet et al.| 
120091 ') . These campaigns led to more than twice the num- 
ber of modes known previously in those two stars, among 
which there were one or two high-order g-modes, and re- 
vealed shortcomings in the excitation predictions for some 
modes. 

A similar excitation problem occurred for 7 Peg, a class 
member found to pulsate in several high-order g-modes 
and low-order p- modes from MOST space-based photome- 
try (H andler et al.l l2009l : Walczak fc Das zyhska-Daszkiewicz| 
bold ). Zdravkov & Pamyatnykh| (|2009l ) suggested that an 
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increase in the opacities by 20% might solve the excitation 
problem of 7 Peg. Another excitation problem was pointed 
ou t for o ne mode of the C0R0T target V1449 Aql (Aerts et 
al.l l201ll ). The recent detection and interpretation of /3 Cep- 
type modes in the 09V pulsator HD 46202 (Briquet et al.| 
1201 ll ) with C0R0T observations showed that the excitation 
problem of massive pulsators is acute, since none of the de- 
tected modes is predicted to be excited for appropriate stel- 
lar models representing this star. In all these studies, the 
effects of rotation on the excitation and amplitude of modes 
were not or only partly treated. Non-adiabatic computations 
taking rotation int o acco unt of the excita tion of modes are 
available (e.g. Leel ll998l ; Townsendl l2005l ) but those of the 
ampl itude of modes are just starting to become available 
fLeel l2012l ). These calculations will be very useful to check 
whether a better agreement between observed and theoret- 
ically excited modes may be found. 

An important result of these asteroseismic studies is 
the evidence for non-rigid interior rotation in some p Cep 
stars. It was shown that V836 Cen, v Eri and 12 Lac rotate 
more rapidly in their inner parts than at their surface. These 
works also allow us to test whether core overshooting has to 
be included in our stellar models. Overshooting represents 
here the amo unt of non-standard mixing processes (Maeder| 
120091 : Mathis| [2uT?j| ). For most of the P Cep targets modelled 
so far, the conclusion was the need for core overshooting 
for a better agreement with the pulsational characteristics. 
The derived core overshooting parameter values are however 
small (around 0.10 local pressure scale heights for V836 Cen 
and HD 46202, and around 0.20 for S Ceti, P CMa, and 
12 Lac) with the exception of 9 Oph for which the value was 
found to be around 0.40. For two cases, v Eri and V1449 
Aql, the core overshooting parameter could be kept to zero. 

Space-based observations revealed other pulsators than 
7 Peg with both g- and p-mode pulsations of SPB and P Cep 
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types (Balona et al. l201ll ). For the CoRoT hybrid HD 50230, 
the observation of deviations from a uniform period spacing 
led to the first exploration of the r egions adjacent to the core 
in a massive star (Degroote et al.l l2010t ). HD 43317 was dis- 
covered by CoRoT to be another hybrid star with a wealth 
of pulsational constraints but also showing the presence of 
chem ical inhomogeneities at its stellar surface (Papics et al.| 
2012). Before space-based asteroseismology, the simultane- 
ous presence of both pulsation and rotational modulation 
was found from ground-b ased s pectroscopy only in j3 Cep 
(Telt ing, Aerts & Mathiasl ll997l ) . V205 2 Oph (Neiner et al.| 
I2003T ), and k Sco (Uytterhoeven et al.l 120051 ). Another im- 
portant discovery by means of space-based CoRoT data is 
the first observation of stochastically excited gravito-inertial 
modes in a massive star. Such modes have recently been de- 
tected in the hot (O9.5-B0) pulsator HD 51452 (Neiner et 
al., submitted), and this type of excitation should probably 
be more widely considered for massive pulsators. 

Among P Cep stars with an asteroseismic modelling, a 
magn etic fie ld has been detected so far in /3 Cep (Henrichs 
et al.l l2000l ; S hibaha shi & Aerts 
(Hubrig et al.l l201ll ; Aerts et al 



2000) and in V1449 Aql 



20111 ). although t he fiel d 



of the latter is still a matter of debate (Shultz et al.l l2012f ). 
In order to understand the effect of a magnetic field on the 
seismic behaviour of /3 Cep stars, a multisite photometric 
and spectroscopic campaign was set up for the equatorial 
star V2052 Oph (my=5.8, spectral type B2IV/V), which is 
known to be a magnetic pulsator with a dominant radial 
mode of frequency 7.14 5 d" 1 and a rotation frequency of 
0.275 d" 1 (Neiner et al.l l2003l ). The star is sligh tly en riched 
in He, revealed a mild N excess (Morel et al.l l2006l ), and 
was considered to be suitable as a seismic target, in view of 
the identified dominant radial mode and the discovery of an 
additi onal low-amplitude non-radial mode by Neiner et al.| 
(|2003f ). 

The multisite photometry of V2 052 O ph is reported in a 
companion paper by Handler et al.| (|2012l ). Here, we present 
the spectroscopic part of the campaign. After a description 
of our dataset in Sect .[21 we discuss the atmospheric parame- 
ters and chemical composition of the star (Sect.[3}. Sect.Uis 
devoted to our frequency analysis and mode identification. 
Besides our line-profile study, we make a comparison with 
stellar models (Sect.0. We end with a discussion about the 
inhibition of mixing by magnetism and conclusions in Sect.[6] 
and Sect .0 respectively. 



2 THE DATA 

A total of 1354 high-resolution spectroscopic exposures was 
assembled using 10 different telescopes spread over both 
hemispheres. In addition, we added the 205 spectra of Neiner 
ct al. (2003, 2012), which allowed us to increase the time 
span of the data set and, thus, to achieve a better frequency 
accuracy. Table[T]summarizes the characteristics of the spec- 
troscopic data. The resolving power A/ A A of the spectro- 
graphs ranged from 20 000 to 70 000. When the instrument 
was not an echelle spectrograph, the spectral domain was 
chosen in order to cover the Si ill triplet around 4567 A be- 
cause these lines are strong without being much affected 
by blending. Moreover, they are dominated by temperature 
broadening so that the intrinsic profile can be modelled with 



a Gaussian, which simplifies the modelling of the line-profile 
variations for mode identification purposes. 

All exposures were subjected to the usual reduction 
process, i.e., we applied debiasing, background subtraction, 
flat-fielding and wavelength calibration. Subsequently, the 
barycentric corrections were determined and the spectra 
were normalized to the continuum by fitting a cubic spline 
function. Special care was taken with regard to the normal- 
isation as it is not only important for an abundance anal- 
ysis but also for a li ne-pr ofile variation study. Indeed, as 
pointed out in Zimsl (|2006l ). a mode identification with the 
FPF method (see Sect. l4.2| l is very sensitive to it. 

The average radial velocity (RV) per site is very similar 
for all spectra, so that there is no evidence that V2052 Oph 
is a spectroscopic binary. To correct for the slightly different 
zero points of the different telescopes, we proceeded in two 
different ways. In a first step, we shifted the spectra in such 
a way that the RV constant of a least-squares sine fit using 
the two dominant frequencies is put to the same value for 
each observatory. In another step, we shifted the spectra so 
that the minimum of the amplitude across the line profile 
for the radial mode (see Fig. [3} is put at the same value for 
each dataset. The two approaches gave the same results in 
our further analysis. 



ATMOSPHERIC PARAMETERS AND 
CHEMICAL COMPOSITION 



Wolff & Heasley| (|l985l ) estimated T eff =23 000 K and 
log g=4.2 based on Stro mgren photometry and fitting of the 
H7 line using Kurucz| \ 19791) atmospher ic mo dels, respec- 
tively. As shown by Nieva & Przybilla| (|2007t ). fitting the 
wings of the Balmer lines using LTE models systematically 
leads to an overestimation of the surface gravit y. On the 
other hand, Niemczura & Daszyhska-Daszkiewicz| (|2005l ) ob- 
tained T e fi=23 350±650 K, log 3=3. 89 and a metallicity of 
-0.25±0.16 dex with respect to solar based on LTE fitting 
of IUE spectr a. Using non-LTE TLUSTY models (Hubeny 
fc Lan z| |l995l ) with a solar helium abundance, Neiner et al.| 
|2003f ) found that T off =21990 K and log #=3. 98 provided 
the best match to three He I lines, but that the Si lines were 
poorly fit in that case. Relaxing the constraint on the He 
content, they inferred T eff =25 200L1100 K, log p=4.20±0.11 
and about a factor two helium overabundance with respect 
to solar (helium abundance by number, y~0.16). A satisfac- 
tory fit to the He and Si line profiles, as well as the contin- 
uum UV flux, was achie ved in that case. 

Morel et al.| (|2006l ) used the non-LTE lin e-form ation 
codes DET AIL/SURFACE (Butler & Giddings| [l985l ; Gid- 
dings| [l98ll ) to carry out an EW-based analysis of the aver- 
age of our 105 exposures from McDonald Observatory. Us- 
ing a series of time-resolved spectra minimises the fact that 
the temperature varies significan tly alo ng the pulsation cycle 
(A?; ff ~900 K; Morton & Hansen lWi Kubiak & Seggewiss| 
1 19841 ) . Fulfilling ionisation balance of Sin/Si ill and fitting 
the wings of four Balmer lines led to T cff =23 000±1000 K 
and lo gfl=4.0dz0.2. Guided by the results of Neiner et al.| 
|2003f ). Kurucz models with He/H=0.178 were used. How- 
ever, this choice has a negligible impact on the atmospheric 
parameters and abundances (Aloge < 0.05 dex). It is not 
possible to further constrain T e ff from ionisation balance of 
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Table 1. Log of the spectroscopic data used in this paper. The Julian Dates are given in days, AT denotes the time-span expressed in 
days, N is the number of spectra and S/N denotes the average signal-to- noise ratio for each observatory measured around the Sim line at 
4552 A. 



\JUacI V clLrWl y 

l 1 N (XlllVl. U 1 L11U llloUl L11I1U11L, 

It'tiUlUllUll, Wd.VUlt_.llH, LI1 1 eLlltiU 111 /l 1 


r 1 iil fic/*nT\f» 
ItlcaLUUc 


Julian Date 
Begin End 


Data amount and 

quality 
AT N S/N 


V_y UScl vcl la J 


2000-2004 




-2450000 










Pic du Midi, France 
(MUSICOS; 35 000; 4489-6619) 


2.0-m TBL 


1730 


3186 


1457 


161 


165 


CN, HH, VG, AH, AT 


2004 




-2450000 










Pico dos Dias Observatory, Brazil 


1.6-m 


3200 


3201 


2 


20 


290 


AB 


(Esp Coudc; 30 000; 4479-4644) 
















Complejo Astronomico El Leoncito Observatory, Argentina 


2.1-m 


3219 


3221 


3 


20 


103 


OP 


(EBASIM; 40 000; 3826-5759) 
















Bohyunsan Astronomical Observatory, Korea 


1.8-m 


3133 


3161 


29 


37 


186 


DM, BL 


(BOES; 50 000; 3751-9803) 
















Dominion Astrophysical Observatory, Canada 


1.2-m 


3134 


3197 


64 


53 


213 


SV 


(45 000; 4457-4603) 
















McDonald Observatory, USA 


2.7-m 


3195 


3198 


4 


105 


215 


GH, ME 


(Coude; 60 000; 3619-10274) 
















Okayama Astrophysical Observatory, Japan 


1.88-m 


3202 


3228 


27 


119 


144 


EK 


Do ouu, oyyi— ^oioj 
















Thiiringcr Landessternwarte Tautenburg, Germany 


2-m 


3142 


3219 


77 


215 


116 


HL 


I^D / UUU, O / UU— 0<±LU) 
















La Silla Observatory, Chile 


1.2-m Euler 


3072 


3282 


211 


60 


104 


KU, KL, JV 


(CORALIE; 50 000; 3876-6820) 
















Mount John University Observatory, New Zealand 


1.0-m 


3154 


3193 


40 


88 


152 


DW 


(HERCULES; 70 000; 4456-7150) 
















Observatorio Astronomico Nacional 
















at San Pedro Martir, Mexico 


2.1-m 


3205 


3211 


7 


637 


282 


RC, JE 


(Echelle Spectrograph; 20 000; 3781-6893) 
















2007-2010 




-2450000 










Pic du Midi, France 
(NARVAL; 65 000; 3694-10484) 


2.0-m TBL 


4286 


5403 


1118 


44 


389 


YF, CM, OT 


Total 








3674 


1559 







neon o wing to the lack of measurable Nell lines (Morel & 
Butlerl l2008t ). 

Sim6n-Dfaz| l|201fj| ) recently reported problems with the 
modelling of some Si lin es using the non-LTE code FAST- 
WI ND (P uis et al.l l2005t ). The model atom used by Morel et 
al.| (|2006t ) is similar in many respects and one may expect 
similar problems to be encountered in our case for some of 
these lines which have been used to constrain T e g (namely 
Sin A4128, which was the only Sin line used, and Sim 
A4813, 4829). We have therefore decided to redetermine T e g 
using lines that are thought to be properly modelled only 
(i.e., Sin A6371 and Sim A4568, 4575, 5740). As can be seen 
in Table(2] this leads to a much reduced abundance scatter 
and a mean value revised upwards to log e(Si)=7.47±0.31 
dex. More importantly, the Si ionisation balance is also sat- 
isfied in that case, thus supporting the previous estimate 
(T cff ~23 000 K). 

Neiner et al.| (|2003h claimed V2052 Oph to be oxy- 
gen weak and helium strong based on spectral fitting us- 
ing TLUSTY models. However, the O abundance they find 
(log e=8.52±0.11 dex) may be regarded as typ ical of the val- 
ues found for other B stars (see Morell l2009l for a review). 



Table 2. Equivalent widths (EWs) rounded off to the nearest 
mA and non-LTE abundances (on the scale in which loge[H]=12) 
for th e silicon lines properly modelled according to Sim6n-Di'az| 
assuming T o fj=23 000 K, logg=4.0 and a microturbulent 
velocity of 1 kms -1 . 



Transition 


EW [mA] 


log e (Si) [dex] 


Sill A6371.4 


43 


7.47 


Sim A4567.8 


100 


7.38 


Sim A4574.8 


66 


7.39 


Sim A5739.7 


85 


7.65 


Mean Sill lines 
Mean Sim lines 




7.47 
7.47 



a Measured on the mean BOES spectrum, as this line falls at 
the extreme edge of an echelle order in the McDonald data and 
is hence unmcasurablc. 

There is instead evidence for a nitrogen excess. The loga- 
rithmic CNO abundances ratios are [N/C]=-0.37±0.08 and 
[N/O]=-0.50±0.13 dex, and appear significantly higher than 



Multisite spectroscopic seismic study of the (3 Cep star V2052 Oph: inhibition of mixing by its magnetic field 5 



the solar values ([N /C]=- O.60±0.08 and [N/O]=-0.86±0.08 
dex; Asplund et al.ll2009l) . This is supported by the results 



of Morel et al.| (|200d ) who found [N/C]=-0.22±0.19 and 
[N/O]=-0.40±0.35 dexQ A nitrogen overabundance is ob- 
served in other magn etic B stars on the m ain seq uence. We 
refer to Morel| (|2010h and Martins et al.| (j2012T ) for a dis- 
cussion on this matter. 

Similarly to Neiner et al.| (|2003j ), a he lium overabun- 
dance was reported by Morel et al.| |2006). However, the 
uncertainties were large and the evidence much weaker in 
that case (?/=0.118±0.032). Although the helium-rich sta- 
tus of V2052 Oph has yet to be unambiguously established, 
this star stands out by the difficulty of fitting the lines using 
a single He abund ance, a fact which explains the large error 
bars (Morel et al.|[200y). This phenomenon might be related 
to surface inhomogeneities/ vertical stratification of helium 
in this magnetic star. The r ecent spectropolarimetric mea- 
surements of Neiner et al.| |2012T l showed that He patches 
are situated close to the magnetic poles. 



4 LINE-PROFILE STUDY 
4.1 Frequency analysis 

For our line-profile study, we considered sufficiently strong 
and unblended lines of different chemical elements (He, C, 
N, O, Mg, Al, Si, S, Fe). These spectral lines were chosen 
among the list of lines used for the abundance analysis study 
by Morel et al. (2006; see their Table A. 1). In particular, we 
analysed the Si III triplet around 4567 A and the He I line 
at 4713 A. The same periodicities were found whatever the 
considered lines. In what follows we only describe our fre- 
quency analysis on the Si ill 4552 A line. 



4-1.1 Equivalent width and radial velocity data 

To perform our frequenc y ana lysis, we used the software 
package FAMIAS0 fZimal l200s1 ). We started by examining 
the radial velocity (RV) and equivalent width (EW) mea- 
surements. In the RV data, we first recovered the known 
pulsation frequency /i =7. 14846 d _1 and its harmonic 2fi 
with an amplitude of 6.7 and 1.5 km s~ , respectively. 
The error estimate of our determined frequency is between 
0.0000 002 d" 1 according to Montgomery & O'Donoghuel 
l|l999l ) (calculated to be a s = y r (6)a s td/ivy r (N)AfAT, 
where <r Bt d is the standard deviation of the final residuals, Aj 
the amplitude of the frequency /, and AT the total timespan 
of the observations) and 0.0003 d" 1 (calculated to be 1/AT). 
In order to be conservative, we adopt five significant digits, 
also for the other frequencies detected in the dataset. The 
EW data also vary with fi. It is due to temperature vari- 
ations caused by the oscillatory displacement of the stellar 



1 These errors are the quadratic sum of the uncertainties in the 
CNO abundances. However, because the abundances have the 
same qualitative dependence against changes in the atmospheric 
parameters, computing the errors on the ratios themselves nar- 
rows down the errors and strengthens the case for an N excess: 
[N/C]=-0.22±0.18 and [N/O]=-0.40±0.22 dex. 

2 FAMIAS has been developed in the framework of the FP6 Euro- 
pean Coordination Action HELAS - |http: //www.hela s-eu. org/ 
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and RV (right) measurements computed from the He I 4713 A 
line (full circles) and the Sin 4552 A line (open circles), after 
prcwhitcning the pulsation frequency /i , for P ro t • 



surface (De Ridder et al. l2002T ). Such a behaviour is observed 
in several /3 Cep stars wit h a hi gh-amplitude mode, such as 
V1449 Aql (Briquet et al.l l2009l and references therein). Be- 
sides fx, we only found a second independent frequency in 
the RV and EW datasets: the frequency 0.27481 d -1 and 
its harmonics. This variability was first discovered from the 
periodic variati on of UV resonance lines observed with IUE 
(Neiner et all 120031 . P la t =3.638833 d) and interpreted as 
due to the stellar rotation. This interpretation is definitely 
confirmed by recent magn etic fie ld data that vary with the 
same period (Neiner et al.l l2012T ). Moreover, as we describe 
below, the observed behaviour in the spectral lines is typical 
of the presence of surface chemical inhomogeneities, which 
leads to line-profile variability as the star rotates. 

In Fig.[T] we show phase diagrams of the EW and RV 
measurements after prewhitening the pulsation frequency 
fx, for Prot- These plots resemble those of the B6 star 
HP 1053 8 2 and of the B7 star HD 131120 (see Briquet et 
al.|[200ll.| [2004Tl . These two stars are He-weak stars whose 
monoperiodic light and line-profile variability is interpreted 
as due to surface chemical spots, which are the result of com- 
plex interactions between radiatively driven dif fusion pro- 
cesses and a stellar magnetic field (Alecian et al.l l201ll ). For 
V2052 Oph, Neiner et al.| (|2003T ) concluded that the ob- 
served rotational modulation is also mainly due to chem- 
ical variations, as the weak magnetic field of the star is 
unlikely to produce strong temperature effects (see discus- 
sion in their Sect. 7.1). As shown in Briquet et al.| (|200ll ). 
a non-sinusoidal RV variation with / ro t and 2/ rot can be 
reproduced by two spots of the same element having a lon- 
gitude difference of ~180° to account for the two maxima 
(minima) separated by a phase difference of ~0.5. More- 
over, the out of phase variation of the EW and RV data 
can be explained if ionized helium is enhanced in regions 
of the stellar surface where Sin is depleted and vice versa. 
This trend was found for both HD 105382 and HD 131120 
by der iving detailed abundance maps for Si and He (Briquet 
et al.l l2004T l. By analogy, the same conclusion may be ob- 
tained for V2052 Oph. An inversion of the observed rotation- 
ally modulated line-profile variations into a two-dimensional 
abundance distri bution (so called Doppler mapping; e.g. 
Piskunov & Ricel ll993l ') could confirm it but is beyond the 
scope of this paper. In the case of V2052 Oph, the difficulty 
is to disentangle the variability due to pulsation and the one 
due to rotational modulation that is also found in spectral 
lines of other chemical elements such as C, N, O and Fe. 
Moreover, while diffusion models in magnetic Ap stars are 
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Figure 2. Amplitude pcriodogram, which is an average of all 
Fourier spectra across the Si III 4552 A profile. The different panels 
show the pcriodograms in successive stages of prew hiteni ng (top 
to bottom) . We follow the notation of Handler et al.| l|2012h for the 
frequencies. Note the different scale of the Y-axis in the different 
panels. 



available for a compa rison with obs erved abundance maps 
(Michaud et al l 198 it Vauclair et al.l ll99ll ). such theoretical 
models need to be developed for stars as hot as V2052 Oph. 



4-1-2 2D analysis 

After a one-dimensional frequency search on the integrated 
quantities, we performed a pixel-by-pixel Fourier analysis 
across the line profile (2D analysis) by means of the tools 
available in FAMIAS. In Fig.0 the mean of all Fourier spec- 
tra across the Si ill 4552 A line is shown for subsequent 
steps of prewhitening. The significance of a frequency cannot 
be derived from the average Fourier spectrum but one has 
to proceed as follows. The wavelength at which the given 
frequency has the highest amplitude is determined. Then, 
the Fourier spectrum at this pixel is computed and the fre- 
quency is retained if its ampli tude exceeds 4 times the signal- 
to-noise (Breger et al.l ll993l ). The noise level is calculated 
in a 5d _1 interval centered on the frequenc y of in terest. 
As for some other P Cep st ars (Telting et all 1 19971 : Schri- 
jvers, Telting & Aertsl l2004l : Briquet et al.l l2005t ). the pixel- 
by-pixel analysis led to additional independent frequencies 
compared to the ID frequency search: /2=7.75603 d _1 and 
/3=6. 82308 d -1 . In the residual periodogram, the frequency 



2/i and combinations of the dominant frequencies /i, f rot 
and 2/ ro t are also seen, as expected in a 2D search (Zima| 
l2008h . Because the 2D analysis is more sensitive to the de- 
tection of high-degree modes than a ID analysis, it indicates 
that /2 and fz are probably not low-degree modes. The four 
frequencies found in our spectra were also d iscove red from 
the multisite photometry by Handler et al.| (|2012t ) and are 
thus certain. 

4.2 Mode identification 

4-2.1 Previous results 

The dominant mode with frequency fi—7. 14846 d _1 was 
first ident ified as radial by Heynderickx, Waelkens & Smey - 
ers| (| 19941 ) and Cugier, Dziembowski & Pamyatnykh| (|l994l ). 
This is unambiguously confirmed by the m ultisit e photomet- 
ric campaign presented in Handler et al.| (|2012T ). As for the 
freque ncies f 2 =7.75603 d" 1 and / 3 =6. 82308 d" 1 , Handler 
et al.| (|2012l ) deduced that both modes correspond to £—i 
or £=6 by comparing the uvy passband amplitudes with the 
theoretically predicted ones. The mo de wi th frequency fz 
was first identified by Neiner et al.| (|2003l ) as an ^3=3 or 
4 mode by me ans o f the MUSICOS spectra. A study by 
Pollard et al.| l|200Sl ) of the HERCULES spectra also re- 
vealed this mode. Their mode identification showed that it 
corresponds to an ^3=4±1 and m,3=2±2 with a best fitting 
obtained for (£$, m 3 )=(4, 2). The combination of these pho- 
tometric and spectroscopic results allows us to safely con- 
clude that 13=4 for f z- Because of geometrical cancellation, 
Handler et al.| (|2012r ) concluded that the most likely iden- 
tification for fa is also ^2=4. We check these conclusions 
through our mode identification aimed at giving constraints 
on (£,m) for /2 and fa (see Sect. l4.2i3l> . 

4-2.2 The method used 

Two independent methods based on spectroscopy are usu- 
ally used to identify the modes of main-sequence pulsators 
hotter t han th e Sun, namely the moment method (Briquet 
& Aert slbOOSt ) and the Fourier parameter fit (FPF) method 
(Zimal 120061 ). Both methods are implemented in the soft- 
ware package FAMIAS and a typical applicatio n of t hem 
on a /3 Cep star can be found in Desmet et al. I |2009l ) for 
12 Lac. For V2052 Oph, /2 and fz are detected in the sec- 
ond moment (line width) but not in the first moment (radial 
velocity). Therefore, we only made use of the FPF method, 
which is appropriate since the projected rotation velocity of 
our target is sufficiently large (t>sini > 20 km s _1 ). 

In the FPF technique, a multi-periodic non-linear least- 
squares fit of sinusoids is computed using all detected fre- 
quencies for every wavelength bin across the profile accord- 
ing to the formula Z + ^ 4 Ai sin[2-K(fit + 4>i)] , where Z is 
the zero-point, and Ai, fa, and (fit are the amplitude, fre- 
quency, and phase of the i-th frequency, respectively. This 
gives, for each pulsation frequency, the observational val- 
ues of the zero-point, amplitude and phase as a function 
of the position in the line profile. The wavenumbers (£, m) 
and other parameters, such as the stellar inclination angle 
i, are then determined in such a way that the theoretically 
computed zero-point, amplitude and phase values derived 
from synthetic line profiles best fit the observed values. The 
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Figure 3. Observed amplitude and phase distributions of the 
/l frequency across the Sim 4552 A line profile (points with 
error bars) compared to the theoretical ones (full line) for 
(l\ , mi)=(0, 0). The amplitude is expressed in units of contin- 
uum and the phase in 2n radians. 



Table 3. The best-fit solutions of our mode identification with 
the FPF method for / 2 =7.75603 d" 1 and / 3 =6.82308 d _1 . 
/l=7. 14846 d —1 is adopted as the radial mode, i and v eq are 
the stellar inclination angle and the equatorial rotational veloc- 
ity, respectively. 



ID 


x 2 


((■2,1712) 




i 












(°) 


(km s^ 1 ) 


1 


1.16 


(4,3) 


(4,2) 


59.7±2.1 


75. Oil. 6 


2 


1.23 


(4,3) 


(4,4) 


33.6±1.9 


118.2±5.1 


3 


1.31 


(4,3) 


(4,3) 


81.6±0.7 


65.3±0.5 


1 


1.42 


(4,2) 


(4,2) 


53.9±3.1 


80.2±3.5 


5 


1.57 


(4,2) 


(4, 4) 


46.9±2.4 


89.1±3.4 


6 


1.63 


(4,2) 


(4,3) 


55.3±7.9 


78.0±7.7 


7 


1.81 


(4,4) 


(4,2) 


45.8±2.6 


90.5±3.5 


8 


1.86 


(4,4) 


(4,3) 


28.5±3.3 


136.1±14.9 


9 


2.01 


(4,4) 


(4,4) 


27.5±1.1 


141.2±5.5 



good ness o f the fit is expressed as a reduced X^-vahie (see 
Zim al l2006h . The approximations made in FAMIAS to com- 
pute the synthetic profiles (see Zima l2006l for details) remain 
valid for our study case. As for stellar rotation, the first or- 
der effects of the Coriolis force are taken into account in the 
displacement field. 



4-2.3 Additional constraints 

We applied the FPF method using the FAMIAS software to 
identify the modes detected in our spectroscopic data. In our 
procedure, we always adopted (4,mi)=0 for the dominant 
mode as deduced from photometry and confirmed by our 
spectroscopy. Indeed, the amplitude and phase behaviour 
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Figure 4. Histograms for the inclination and equatorial rota- 
tional velocity of the star derived from the FPF method, for the 
combination (I2, m2)=(4, 3) and (£3, m 3 )=(4, 2). 



across the line profile for this mode is very typical for a 
radial pulsation mode (see Fig.[3jl . 

We adopted the values M = 9 Mq and R — 5 Rp> fo r 
the mass and radius, respectively (Handler et al.l 120121 ). 
The mode identification results are robust when using dif- 
ferent values for these parameters, within the errors. Be- 
fore starting the mode identification, the parameter space 
was restricted by estimating values for the projected ro- 
tational velocity wsini, the width of the intrinsic gaus- 
sian profile cr, and the equivalent width EW from a least- 
squares fit of a rotationally broadened synthetic profile to 
the zero-point profile. A genetic optimization with the fol- 
lowing free parameters was afterwards adopted: £ £ [0, 8] 
with a step of 1, m € [— I, £] with a step of 1, the sur- 
face velocity amplitude a £ [l,50]kms _1 with a step of 
lkms -1 , the stellar inclination angle i £ [1,90]° with a 
step of 1 °, wsini G [50, 70] kms -1 with a step of 0.1 kms -1 
and a € [5, 15] kms -1 with a step of O.lkms -1 . The fitting 
between the observed and theoretical zero-point, amplitude 
and phase was carried out by applying genetic optimization 
routines. First, we used mono-mode fits in order to speed up 
the computations. Once constraints on the parameters are 
obtained, we performed a multi-mode fit. 

The results clearly exclude £ below 3 and above 5 for 
both low-amplitude modes. Combining this out come with 
the photometric one (£=4 or 6; Handler et al.l I2012T I. we 
find that £=4 for / 2 and fs . In order to fit the phase across 
the profile, we also found that m > (where a positive m- 
value denotes a prograde mode) for both modes. Moreover, 
a visual inspection of the solutions revealed that we could 
also exclude the (£, m)-combinations involving m=l. The 
matching between the observed and synthetic amplitudes 
across the profile was not satisfactory in these cases. 

In Tabled we list the remaining nine acceptable (£, re- 
combinations and give estimates of the stellar inclination an- 
gle i and equatorial rotational velocity u eq obtained through 
a simultaneous fit of the modes. Histograms, constructed 
using the x 2-va l ue as a weight, for the solution with the 
lowest x 2 ( see Table|3} are displayed in Fig.[4] Fig. [5] de- 
picts the matching between the observed and theoretical 
amplitude and phase across the line, for the combination 
(&,m 2 )=(4,3) and {£ 3 , m 3 ) = (4, 2). 



5 COMPARISON WITH STELLAR MODELS 

We computed stellar models accounting for the derived ob- 
servational properties of V2052 Oph. To this end, we used 
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f 2 =7.75603 d" 



f 3 =6.82308 d" 




Velocity (km s 



Figure 5. Amplitude and phase distributions for fa and fo for the Sim 4552 A line (points with error bars) and comparison with the 
theoretical ones (full lines) for the solution with (£2 , rri2)=(4, 3) and (£3, rri3)=(4, 2). The amplitudes arc expressed in units of continuum 
and the phases in 27r radians. 



the evolutionary code CLES (Code Liegeois d'Evolution 
Stellaire; Scufiaire et al.l l2008al ) with the input physics de- 
scribed in Briquet et al.| (|201ll ). For each stellar model con- 
sidered, the theoretical pulsation frequency spectrum was 
calculated in the adiab atic ap proximation with the code 
LOSC (Scufiaire et al.l l2008bl ). Once models that fit the 
observed modes were selected, we checked the excitation 
of the pulsatio n mod es with the lin ear no n-adiabatic code 
MAD (Dupretl l200ll ; Dupret et al.l 120021 1. In our compu- 
tations, the effects of rotation on the pulsation frequencies 
were taken into account to first order, which is sufficient 
for our approach as justified a posteriori at the end of this 
section. 

By varying the mass (M G [7.6, 20.0] Mq with a step 
of 0.1 M Q ), the hydrogen mass fraction (A' G [0.68, 0.74] 
with a step of 0.02), the metallicity (Z G [0.010, 0.018] 
with a step of 0.002) and the core overshooting parameter 
(q v G [0, 0.5] with a step of 0.05 - expressed in local pres- 
sure scale heights), we computed models from the zero-age 
main-sequence (ZAMS) to the terminal-age main-sequence 
(TAMS). 

Given the rotation period of the star and the model 
radius when /1 is identified as the fundamental, we get an 
equatorial rotational velocity between 71 and 75 km s _1 for 
V2052 Oph. With vsini G [61,65] km s"\ the star has an 
inclination i G [54,66]°. This range of inclination angles is 
in full agreement with the one determined by Neiner et al.| 
(|2012l ) who obtained an inclination between 53° and 77° by 
synthesizing the Stokes V profiles of the star with a centered 



or offcentered dipole model. Stellar models with the radial 
mode as the first overtone (or higher overtones) have larger 
radii, implying an inclination lower than 50° that is not 
compatible with the constraints deduce d from the magnetic 
field data. We refer to Hander et al.| (|2012l ) for a deeper 
discussion. 

From the models that we computed, we thus selected 
the ones along the main-sequence evolutionary track that 
fit the radial fundamental mode to within 0.05 d _1 . The po- 
sition of these models in the log Teff-log g diagram is com- 
pared to the spectroscopic error box of the star in Fig. [5] 
For all models matching /1, we computed the theoretical 
frequency spectrum of low-order p- and g-modes with a de- 
gree of £=4. The rotationally split frequencies were derived 
using the Ledoux constant as /„f m = fnto + mftni] 'rot, where 
/3 n e is a structure constant depending on the stellar model. 
Afterwards, to each model fitting the radial fundamental 
mode, we assigned a x 2 -value, which compares the observed 
frequencies /° bs with the theoretically computed ones /* h , 
as follows: 

3 /robs rth\2 

where a is taken to be 0.05 d _1 to account for uncertain- 
ties in our modelling prescriptions (see last paragraph of 
this section). The models corresponding to one of the nine 
(£, m)-combinations of Table|3] and having a X/ < 1 were 
retained, also requiring the excitation of the radial mode. 
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Table 4. Physical parameters of the models fitting the three pulsation modes detected in V2052 Oph according to their mode identifi- 
cations listed in Tablc[3] using the same IDs. The last line corresponds to the parameter range of our best models for the star, i.e., with 
a Xf as defined in Eq. JTJ lower than 1 and excluding solutions with ID 7, 8 and 9. 



ID 


mass 


radius 


Z 






log g 


log(L/L ) 


age 






(M Q ) 


(R©) 






(K) 






(Myr) 




1 


[8.2, 9.0] 


[5.16, 5.31] 


[0.012, 0.016] 


[0.00, 0.15] 


[21370, 22190] 


[3.93, 3.94] 


[3.72, 3.76] 


[19.9, 23.7] 


[0.28, 0.32] 


2 


no model 


















3 


[8.9, 9.6] 


[5.26, 5.38] 


[0.010, 0.012] 


[0.00, 0.10] 


[22790, 23570] 


[3.95, 3.96] 


[3.85, 3.89] 


[17.3, 18.3] 


[0.27, 0.31] 


1 


[8.3, 9.0] 


[5.17, 5.30] 


[0.012, 0.016] 


[0.00, 0.10] 


[21600, 22250] 


[3.93, 3.94] 


[3.73, 3.77] 


[19.2, 21.9] 


[0.27, 0.30] 


5 


no model 


















6 


[9.0, 9.1] 


[5.27, 5.29] 


0.010 


0.05 


[23320, 23620] 


3.95 


[3.89, 3.87] 


[16.9, 18.0] 


[0.25, 0.27] 


7 


[8.1, 9.0] 


[5.15, 5.31] 


[0.016, 0.018] 


[0.00, 0.15] 


[20550, 21720] 


[3.92, 3.94] 


[3.64, 3.74] 


[18.7, 26.5] 


[0.30, 0.35] 


8 


[7.8, 9.6] 


[5.10, 5.39] 


[0.010, 0.018] 


[0.00, 0.35] 


[20260, 23680] 


[3.91, 3.96] 


[3.61, 3.90] 


[15.6, 28.9] 


[0.28, 0.37] 


9 


[7.7, 8.2] 


[5.10, 5.19] 


[0.014, 0.018] 


[0.25, 0.40] 


[20350, 21620] 


[3.91, 3.92] 


[3.62, 3.71] 


[24.5, 30.1] 


[0.33, 0.38] 




[8.2, 9.6] 


[5.16, 5.38] 


[0.010, 0.016] 


[0.00, 0.15] 


[21370, 23620] 


[3.93, 3.96] 


[3.72, 3.89] 


[16.9, 23.7] 


[0.25, 0.32] 



These models are displayed as crosses in Fig.[6] and their 
physical parameters are listed in Table|3] 

Although we permitted the two £=4 modes to have the 
same radial overtone, none of the models matching our ob- 
servations are such that and /3 belong to the same ro- 
tationally split multiplet. In all models appropriate for the 
star, the frequencies /2 and f$ have ni=— 3 (.gs-mode) and 
713 = — 2 ((72-mode), respectively. Both frequencies are pre- 
dicted to be excited by non-adiabatic computations in all 
cases. In fact, a scan of the modes with £=4 shows that only 
those with n=— 1, —2, and —3 are theoretically excited. Fol- 
lowing theory, modes with n=—l, 1 for £=1, n=— 1, —2 for 
£—2, and n=—l, —2, —3 for £=3 are also excited in the same 
frequency range as the observed modes. Clearly, if they are 
present in V2052 Oph, their amplitudes are lower than the 
observed £=4 modes. It may be explained by the rotation 
rate of the star that is somewhat higher than in other well- 
studied ft Ce p star s with dominant low-degree modes (see 
Hander et al.l l2012l for more comments). 

In Table[3] the inclination angles i of solutions with ID 
1, 4 and 6 are fully compatible with the stellar modelling. 
Solutions with ID 2 and 5 are rejected because no model 
with the corresponding (£, m)-combinations could be found 
(see TableQ. If one limits the inclination range around 60° 
during the mode identification as performed in Sect. 14.21 so- 
lutions with ID 7, 8 and 9 have even higher x 2 -values so that 
we can exclude them. However, solution with ID 3 is kept 
since the x 2 -vahie is similar as before (1.33). The most prob- 
able identification for both /a and is thus (£, m)=(4, 2) or 
(4,3). 

The parameter ranges of our best models for V2052 Oph 
are listed in the last line of Tabled They are obtained by 
only considering the solutions with ID 1, 3, 4 and 6 and 
are represented as white circles in Fig. [6] In this figure, we 
also distinguish the models associated to ID 1 as black cir- 
cles. An important result is that only stellar models with 
no or mild core overshooting can explain our observational 
constraints. This is further discussed in the following sec- 
tion. One also notices the excellent agreement between our 
best fitting model parameters and the T e g and log g- values 
derived from spectroscopy. Among ft Cep stars with an as- 
teroseismic modelling, suc h a co nsistency was also obtained 
for 12 Lac (Desmet et al.l l2009l ) but it is not often the case 



(e.g. Briquet et al.l l2007l ; Briquet et al.l l201ll ; Aerts et al] 

Eon]). 

In order to assess the influence of rotation on the pulsa- 
tion frequency values, we calculated, usin g the 2D method 
described in Reese, Lignieres & Rieutord| (|2006r i, the funda- 
mental radial mode and the £—4, n——2, —3 multiplets of ro- 
tating polytropic models with a polytropic index of 3 (which 
is typical of a radiative zone), and a mass and polar radius 
set to 9M0 and 5R©, respectively. A comparison between a 
first order approximation of the frequencies and the 2D cal- 
culations yielded differences around 0.05 d _1 for the radial 
mode, and 0.02 d _1 or less for the £=4 modes, for an equato- 
rial ro tation al velocity of 75kms~ 1 . As explained in Ballot 
et al.| (|2010l ). p- modes are more affected by the centrifugal 
deformation, given that their energy is concentrated towards 
the outer layers of the star, thereby explaining the larger 
difference for the radial fundamental mode. Although these 
differences can be detected, given the length of our obser- 
vational run, they remain below the percent level compared 
to the frequencies. Furthermore, they are much smaller than 
the distance between the modes (n, £, m) = (—2,4,0) and 
(—3, 4, 0), and the spacing between consecutive members of 
the g-mode multiplets when the frequencies are expressed in 
an inertial frame. The above differences explain why we con- 
sidered an error bar of 0.05 d _1 on the pulsation frequencies 
for our comparisons between observed and model frequen- 
cies. 



6 DISCUSSION 

V2052 Oph is very similar in terms of stellar parameters to 
the ft Cep star 9 Oph, except that V2052 Oph has a higher 
surface rotation velocity (see Table[5|. 9 Oph rotates at 6% 
of its critical breakup velocity, while V2052 Oph rotates at 
16%. For such stars, we expect a larger core overshooting 
param eter for the faster rotator. Indeed, Browning et al.| 
l|2004l ) showed that when the rotation rate is increased (for 
given stratification and diffusivities) the penetration of flows 
coming from the convective core in the surrounding stably- 
stratified radiative envelope becomes more intense because 
of the stabilisation of the flow by the Coriolis acceleration. 
In addition, the rotational mixing associated to the possible 
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Figure 6. The error box represents the position of V2052 Oph, 
deduced from spectroscopy, in the log T c ff-log g diagram. The 1, 2 
and 3-cr regions in T c g are shown. The light grey band corresponds 
to models that fit the radial fundamental mode. Our best models 
for the star are displayed as crosses (considering all nine (£, re- 
combinations of Tablcl4t . white circles (with all solutions but ID 
7, 8 and 9), black circles (for the solution with ID 1). We refer to 
the text for more explanations. 



differ ential rotation increa ses with rotation rate (Mathis & 
Zahnl l2004l ; Maederl [20091 ). 

The reverse result is obtained since the core overshoot- 
ing parameter is lar ger fo r Oph than for V2052 Oph. In- 
deed, Briquet et al.| (|2007h obtained « OV =0.44±0.07 H p with 
Z=0.009 to 0.015 for 9 Oph, using the same tools and input 
physics as here. Using a different stellar evolution modelling 
code an d a 2D linear adiabatic pulsation code, Lovekin & 
Goupil| (|2010h found a ov =0.28±0.05 H p assuming Z=0.02. 
The results of both studies are compatibl e since increasing Z 
decreases a ov , as shown in Briquet et al. Such values 

were also fo und in the modelling by Daszyhska-Daszkiewicz 
& Walczak| (|2009m . 

The inhibition of mixing as derived for V2052 Oph can 
be caused by processes that are able to damp differential ro- 
tation in the stellar radiation zone, i.e. either st ochast ically 
excited internal waves (e.g. Talon fc Charbonnell l2005l ; Pan- 
tillon Talon & Charbonnell l2007T ) or a magnetic field (e.g. 
Mossl ll992T l. The relative size of the convective core of the 
two stars is similar: 12.5% of the total radius for V2052 Oph 
and 13.5% for Oph. Therefore, if internal waves transport- 
ing angular momentum are excited by the turbulent con- 
vection in the core, they are expected to have the same net 
effect in the two stars. V2052 Oph is known to be a mag- 
netic star ( Neiner et al.l l2003h while 9 Oph is not (Hubrig 
et all 120061 '). V2052 Oph hosts a fossil m agnet ic field with 
Bpoi of the order of 400 G (Neiner et al.l l2012f ). Therefore, 
the presence of the magnetic field in V2052 Oph seems to 
be the most likely explanation for the low core overshoot- 
ing parameter, which represents the non-standard mixing 
processes. We test below whether the observed field can be 
sufficient to inhibit mixing in V2052 Oph. 

W e first used the criterion from Eq. (3.6) in Zahn| 
(|201ll ), which is based on the equation of angular momen- 
tum transport taking into account the magnetic field derived 



Table 5. Comparison of the parameters of V2052 Oph and 6 Oph. 



Parameter 


e Oph 


V2052 Oph 


T eff (K) 


22260±280 


22500±1100 


logg (dex) 


3.95±0.01 


3.95±0.02 


M (M Q ) 


8.2±0.3 


8.9±0.7 


x c 


0.38±0.02 


0.29±0.04 


z 


[0.009, 0.015] 


[0.010, 0.016] 




0.44±0.07 


0.07±0.08 


V e q (km s™ 1 ) 


29±7 


73±2 


f2/f2 cr it 


0.09±0.03 


0.23±0.01 



The parameters for 6 Oph are taken from Briquet et al.| l l2007h . 

by Mathis & Zahn| (|2005l ): 

Bl it =^pp±. (2) 

Jam 

In this equation, we used the time already spent by 
V2052 Oph on the main sequence, its radius and mass, as 
well as the surface rotation velocity. We found that the mean 
critical field strength in the radiative zone for suppressing 
differential rotation in this star is B C rit~70 G. According 
to the ratio (30) between internal and surface fields derived 
by Braithwaite (2008; see Fig. 8 therein) , this corresponds 
to a critical field strength at the surface of -B C rit,surf=2 G. 
Although the Zahn criterion is an approximation and we 
took global parameter values in our numerical application, 
it provides a good idea that a very weak fossil magnetic 
field is enough to inhibit possible initial differential rotation 
and thus rotational mixing in V2052 Oph or similar early 
B stars. Thus the surface magnetic field of 400 G present 
in V2052 Oph is strong enough and V2052 Oph is rotat- 
ing uniformly considering the approximations of the Zahn 
criterion. 

We then co nsider ed the criterion from Eq. (22) de- 
fined by Spruit| (|1999T ), which provides the critical initial 
field strength above which the magnetic field remains non- 
axisymmetric and rotation becomes uniform: 

In this equation we consider the radiative zone. From our 
CLES models we derived that this radiative zone encom- 
passes 81% of the stellar mass and its average density is 
0.07 g cm" 3 . In addition, we assumed that the initial rota- 
tional velocity in the radiative zone is a few times the current 
surface rotation velocity. We used a magnetic diffusivity of 
2 x 10 6 , typical of the radiative zone of a B2 star (Brun, pri- 
vate communication; see also Augustson, Brun & Toomre| 
l201ll ) and a degr ee of differential rotation (q = r\ VT2|/fl ; 
Eq. (13) of Spruitl ll999f l equal to 1 following Spruit | (|l999h . 
We obtain that the initial critical field is j3 C rit,i nit^4 30 G. 
Using again the ratio provided by Braithwaite| (|2008l ). but 
this time between the radiative zone and the surface, i.e. 
typically a factor 10, the critical initial surface field is 
S cr it,init,surf~40 G. The corresponding current critical sur- 
face field is below this value, because of dissipation caused 
by the ohmic diffusion and the uniform rotation state caused 
by the Lorentz torque. Again, the Spruit criterion is based 
on approximations. Nevertheless, it shows that the surface 
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magnetic field observed in V2052 Oph (400 G) is sufficient 
to explain the observed non-axisymmetric magnetic config- 
uration and lack of mixing associated to uniform rotation. 

More precise criteria extracted from detailed modelling, 
in particular taking into account the full coupled transport 
processes and the geome trical configuration (for example fol- 
lowing Mathis & Zahn l [200a ; Spada, Lanzafame & Lanzal 
20ld : Duez & Mathisl 120101 Duez, Braithwaite & Mathis| 



20101 ). are required to determine more precisely the level of 



magnetic field necessary to inhibit mixing in B stars. Never- 
theless, our first estimate with both criteria above confirm 
that it is probably the presence of a fossil magnetic field in 
V2052 Oph that reduces the core overshooting parameter, 
i.e. the amount of non-standard mixing processes. 

As a consequence of our findings, the usual interpre- 
tation that rotational mixing is responsible for the N- 
enric hment observed in early B-type stars (Meynet et al.| 
l201ll ) cannot be applie d to V2052 Oph. Preliminary results 
of Bourge et al. (2007) showed that radiatively-driven mi- 
croscopic diffusion may be an alternative explanation for 
this abundance peculiarity in V2052 Oph and other targets. 
Complete diffusion computations including radiative forces 
in /3 Cep-type stars, which are not yet available, should in- 
clude magnetic fields as there seems to be a higher incidence 
of nitrogen overabundance in stars with a detected magnetic 
field of the same order (i.e. a few hun dred G auss) as the one 
detected in V2052 Oph (Morel et al.l l2008h . 



7 CONCLUSIONS 

Our study was based on extensive ground-based, high- 
resolution, high-S/N, multisite spectroscopic measurements 
spread over many years. Similar efforts were already per- 
formed for two other ft Cep stars, u Eri and 12 Lac, lead- 
ing to around ten independent pulsation frequencies. For 
V2052 Oph, our campaign only led to the detection of three 
pulsation modes but, contrary to the two other cases, the 
rotation period is also found, which provides a strong ad- 
ditional constraint. The difficulty to detect more modes in 
V2052 Oph may be explained by lower amplitudes of the 
modes due to a higher rotation rate compared to the two 
other stars. 

Comparing our identifications of the degrees of the 
modes with the ones coming from the photometry, we are 
confident that the dominant mode is radial while the two 
additional modes have t=A. Moreover, combining the out- 
come of the FPF method with that of our basic mod- 
elling, we conclude that the non-radial modes are prograde 
(£, m)=(4,2) or (4,3) components of two multiplets with 
n=— 2 and n=— 3. The surface equatorial rotational velocity 
was also deduced 73±2 km s . 

The observation of three independent modes with 
constraints on their wavenumbers along with the rota- 
tion period was enough to deduce ranges for the mass 
(M e [8.2,9.6] Mq), central hydrogen abundance (X c 6 
[0.25,0.32]) and other parameters of V2052 Oph that are 
furthermore fully compatible with the spectroscopic T e g 
and logy, and with the inclination angle i determined from 
the modelling of the magnetic field. More importantly, we 
showed that only models with no or mild core overshooting 
could satisfactorily reproduce our observational properties. 



We made a comparison with 6 Oph because it is an 
asteroseismically modelled ft Cep star having similar funda- 
mental parameters as our object of interest. We would ex- 
pect a larger core overshooting parameter value in the faster 
rotator but the contrary is found. It can be explained by the 
fact that V2052 Oph is a magnetic star while 6 Oph is not. 
Indeed, using two approximate criteria, we showed that the 
magnetic field present in V2052 Oph is strong enough to in- 
hibit non-standard mixing processes in its interior. In order 
to check if the internal rotation of V2052 Oph is indeed uni- 
form or not, we would need additional constraints such as 
those obtained from a high-precision light curve as given by 
the C0R0T and Kepler satellites. Indeed, space-based pho- 
tometry has already revealed many pulsation frequencies, 
even in objects that wer e thou ght to be mono period ic from 
the ground (Aerts et al.l l2006l ; Degroote et al.l l2009l ). 
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